Mon. Not. R. Astron. Soc. 000, [|-| (2001) Printed 1 February 2008 (MN WI^ style file v2.2) 



A practical theorem on gravitational wave backgrounds 
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1 INTRODUCTION 



ABSTRACT 

There is an extremely simple relationship between the spectrum of the gravitational 
wave background produced by a cosmological distribution of discrete gravitational 
wave sources, the total time-integrated energy spectrum of an individual source, and 
the present-day comoving number density of remnants. Stated in this way, the back- 
ground is entirely independent of the cosmology, and only weakly dependent on the 
evolutionary history of the sources. This relationship allows one easily to compute the 
amplitude and spectrum of cosmic gravitational wave backgrounds from a broad range 
of astrophysical sources, and to evaluate the uncertainties therein. 

Key words: gravitational waves - diffuse radiation ~ binaries: close - black hole 
physics - relativity 



The strongest discrete sources of gravitational waves are those which radiate large amounts of energy in a time (short 
compared to the age of the universe) before or after a catastrophic event. Examples include supernovae, accretion-induced 
collapse, black-hole formation events and merging compact binary systems, which may involve white dwarfs, neutron stars or 
black holes of the stellar or super-massive sort. It is becoming increasingly important to understand the 'brightness' of the 
night sky in gravitational waves due to the cosmic superposition of such sources. This is crucial to the design of LISA, the 
NASA/ESA Laser Interferometer Space Antenna, and for proposed follow-on missions at higher and lower frequencies which 
might search for primordial stochastic backgrounds, e.g. from inflation or phase transitions in the early universe. There have 



spectra or waveforms of sources with assumed evolutionary histories over cosmological volu mes (Coward, Burman & Blair 



2001 



Schneider et al 2001 



Ferrari, Matarrese & Schneider 1999 



Rajagopal & Romani 1995) 



We show here that there is an extremely simple relationship between the spectrum of the gravitational wave background 
produced by a cosmological distribution of such sources and the present-day comoving number density of remnants. The 
background is entirely independent of the cosmology, and in most cases is almost independent of the evolutionary history of 
the sources. This simple relation allows one to evaluate quickly uncert ainties in estimates of gravitational wave backgrounds 



and to survey new backgrounds. A companion paper ( Phinney 2001 ) in this way surveys the gravitational wave sky, and 



points out a number of previously unrecognized backgrounds of potential significance for LISA and future missions. 

In section ^ we give a simple physical derivation of the theorem. Section ^applies the theorem to the backgrounds produced 
by non-relativistic binaries in circular orbits, and compares to previous calculations of the cosmic background from double- 
degenerate binaries and merging super-massive black holes in galactic nuclei. The reader wanting quick numbers for other 
circular orbit sources should use equations p^|-p^. Section ^ presents a formal derivation of the theorem, and shows explicitly 
that it is valid for all source lifetimes, including ones much longer and much shorter than the duration of a measurement 
experiment. 



2 THE THEOREM: PHYSICAL DERIVATION 

Let fr be the frequency of gravitational waves in a source's cosmic rest frame, and / the frequency of those waves observed 
today on earth, fr = /(I -|- z). Let the total outgoing energy emitted in gravitational waves between frequency /,■ and fr -\- dfr 
be 
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dEgyj , , 

IT ■ 

This energy, like fr is measured in the source's cosmic rest frame, and is integrated over aU sohd angles and over the entire 
radiating lifetime of the source. 

Let the number of events in unit comoving volume which occur between redshift z and z + dzhe N{z)dz. Define Qg^^if) 
to be the present-day energy density per logarithmic frequency interval, in gravitational waves of frequency /, divided by 
the rest-mass energy density pcC^ that would be required to close the universe. Then the total present day energy density in 
gravitational radiation is 

Pcc'ilgAf)df/f^ 1^°° l^fh'M)^ , (2) 

where pc — 3Hq / (SivG) is the critical density of the universe, and he is the characteristic amplitude of the gravitational wave 
spectrum over a logarithmic frequency interval din/ = df /f . he is related to the one-sided (0 < / < oo) spectral density 



Sh.i of the gravitational wave background (cf. Thorne (1987)) by hj.{f) — fSh,i{f )- If the spectral density 5^,2 is defined 
to be two-sided (-oo < f < oo, cf. UngareUi & Vecchio (pOOl|), Cornish (|200l|)), the relation is hl{f) = 2fSh.2{f)- ^l^'^ is 



often called the strain noise, since the mean square signal strain output of an interferometer is simply the integral over all 
frequencies of Sh times the response function TZ{f) of t he in terferometer (of order unity for wavelengths much longer than 



the interferometer -cf. section 3.2 of Cornish & Larson (2001)). 

In any homogeneous and isotropic universe, the present-day energy density Eg^ must be equal to sum of the energy 
densities radiated at each redshift, divided by (1 -|- z) to account for the redshifting of the gravitons since emission; 

Equating the two expressions in equations (^) and (^) for Egu, frequency by frequency, we find 



1 f . dEa 



dz . (5) 

This theorem is our principal result. It has the simple physical interpretation that the energy density in gravitational waves 
per log frequency interval is equal to the comoving number density of event remnants, times the (redshifted) energy each 
event produced per log frequency interval. Notice that the theorem does not depend upon the cosmological model, except 
for the assumption of a homogeneous and isotropic universe. Nor does it depend on any property of the sources (beaming, 
polarization, etc) except for their time-integrated energy spectrum, provided they are randomly oriented with respect to earth. 
If more than one type of source is important, the right hand side of equation ^ should be summed over the source types i, i.e. 

dEgW \ . 1 J K \ ^ AT I \ ( f dEgtij^i 



N[z) yfr^j is replaced by }_^N,{z) (^,A~^j • (6) 
The a lert re ader will detect some kinship between the frequency-integrated version of this theorem (equation W) and that 



of Soltan (1982) relating the contribution of quasars to the electromagnetic brightness of the night sky to the local space 
density of remnant super-massive black holes. This relation becomes even closer if the source we are considering is black 
holes or other compact objects of mass M2 gravitationally captured in circular orbits by super-massive black holes of mass 
Ml ^ M2. For them J {dEgu,/dfr) dfr = M2C?eb, where eh = {eisco ~ / dt), where dsco is the dimensionless binding energy 
per unit rest mass of the innermost stable circular orbit (0.057 for Schwarzschild black holes, 0.42 and 0.038 respectively for 
prograde and retrograde equatorial orbits about maximally rotating Kerr black holes) , and J ch dt is the (small) portion of 
the binding energy radiated down the black hole horizon (note that this can be, and is, negative for prograde orbits around 
rapidly rotating black holes: the orbiting mass extracts some of their rotational energy -see table VII of Finn & Thorne 



(2000)). The change in the mass of the capturing black hole is AA/i = h'hil — ej). Thus 



Eg^^ iV(^)-^-^AA/lC^d^ = p.c^/™(-^-^), (7) 

where fm is the fraction of the present-day comoving mass density p, in super-massive black holes which was grown by the 
capture of compact objects (as opposed, say to growth by axisymmetric accretion of gas which results in negligible gravitational 
radiation) . 

In section H we give a more mathematical proof of the theorem, but we first give an example of its use. 



3 AN EXAMPLE: MERGING BINARIES IN CIRCULAR ORBITS 

In the Newtonian limit, a circular binary of component masses Mi and M2 which merges due to gravitational radiation losses 
in less than the age of the universe has 
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dE,^ _ n 1 {GMf^' farf<f<f (H) 

where A4 is the chirp mass, A4''^^ = M\M2{M\ + M2)~^^'^ ■ Notice that to derive equation one need know nothing about 
gravitational radiation except that it has twice the orbital frequency, and that the orbital binding energy is removed by the 
gravitational radiation. 

The lower limit /min is set by the separation of the system at its birth (or circularisation, whichever comes first), and 
our derivation applies only to those systems whose initial separation is small enough that their time to merge is much shorter 
than the Hubble time at their birth. For example, a pair of 0.3Mq white dwarfs must have /min > 10~* Hz in order to merge 
through gravitational radiation in < 10^° y. Over their lifetime, merging systems radiate a broad spectrum of frequencies, 
up to an upper limit /max set by the frequency at which the two bodies come into Roche lobe contact, or at which tidal 
dissipation begins to dominate gravitational radiation in determining the orbital evolution. For a binary whose least massive 
s tar is a white dwarf of mass M2 = M^dMQ, one can use the fit to white dwarf mass-radii of equation (17) of Tout et al 
(1997) and the approximation (good to better than 10%) that the orbital frequency Qb = 7r/r of a binary is related to the 
equivalent volume radius of M2's Roche lobe Rl2 by = O.IGM2/ Rj^2j to find the maximum gravitational wave frequency 
emitted before the onset of mass transfer: 

- °-°^ ^i - (aWi.Vt/- • 



For pairs of neutron stars, /max — 1.4 x 10 Hz ( Uryu, Shibata fc Eriguchi 200C ; Shibata fc Uryu 2000 ). For small bodies 
spiraling into a non-rotating black hole of mass more detail and the 

generalization of equation pi to fully relativistic orbits around rotating black holes, see Phinney (2001)). Systems whose initial 
separation is so large that their time to merge is much longer than the Hubble time at their birth will have /max — /min, and 
are neglected here. 

Inserting equation |^ into equation ^ gives the gravitational wave background at / < /max from a population of inspiraling 



binaries 

9 c2R 



^,Uf) = ^:^(GMf' f^'No{{l + z)-'^') s.nd (10) 




4 1 {GMfr^ , 



hiif) = ^-2^Jl/^M(l + ^r (11) 
where A^'o = JJ^ N{z) dz is the present-day comoving number density of merged remnants, and 

— max[0, /min// — 1] and ^^max — /max// 1 Can be Set respectively to and oo except for / just below /min or /max. In 
terms of the merger rate per comoving volume N{tr), N{z) — Ndt^/dz, where 

-7- ~ rr — ^ s N ) and (13) 

dz Hoil + z)E{z) ^ ' 

E{z) = {QM(l + zf + nk{l + zf + nAy^^. (14) 



Currently favoured values are D.m = 0.33, = 0.67, Qk = 0, Ho = 65 km s"^Mpc"^ (cf. Netterfield et al ( |200l[ )). The 
value of ((1 + z)~^^^) is not very sensitive to the details of N{z). For example, if we take N increasing rapidly in the past, 
propo rtional to the cosmic star formation rate as a function of redshift given in equation (6) of Madau, Haardt & Pozzetti 
(2001) in a fiat = 1 universe, we have {{1 + z)~^^'^) — 0.74, while a time-independent iV in a fiat Qa = 0.67 universe gives 
{(l-l-z)-i/^) =0.80. 

The one approximation made in equations ^ or |ll| is that the lifetime has been assumed short compared to the expansion 
time of the universe -i.e. all of the inspiral occurs at the same redshift as the merger. We are thus neglecting the differential 
redshift of the very earliest parts of the inspiral which occur at higher redshift than the merger. But because of the tmrg oc a** 
dependence on initial separation a, sources generally either merge quickly or not at all, so this approximation is poor only for 
a small percentage of sources with finely tuned parameters. 

In convenient numerical forms, equations |l^ and ^ and the 1-sided strain noise Sj/'^{f > 0) = ftc/vT (the two-sided 

Sl^2i~°° < / < °o) is Sl^^ /V^) become 



Mq J vio-=^Hzy ^Mpc-^^y 0.74 

Mq) ) \ 0.74 



M/) = 3.0xl0-(^) (,^) (^) ( ^^--^. ' ] , (16) 
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The redshift averages in equations |15|-|17| are frequency-independent constants for /min/(l + z») ^ f ^ /max/(l + z*), where 
is the median source redshift; outside that range the edge effects of equation |l2| become important. 
We now use these to estimate the cosmic gravitational wave background from merging double-degenerate binaries. For 
our Ho = 65 km s~^Mpc~^, Fukugita, Hogan & Peebles (1995) give the mass density of sta rs in the universe as = 
4.4 X IO^MqMpc"^. After 10^°y a stellar population with the initial mass function of Kroupa (2001) has 0.28 white dwarf 
remnants per solar mass of stellar material, giving a present-day comoving number density of 1.2 x lO^WD Mpc"^. Of 
these white dwarfs, under standard assumptions of population synthesis models (equal numbers of binaries per log initial 
orbital separation, binary fraction 0.5, flat binary mass ratio distribution), a fraction 0.015 of the white dwarfs will undergo 
common envelope evolution leading to a pair of ~ 0.3Mq white dwarfs which will merge in < lO^^y. Thus we estimate 
iVo = 1.8 X lO'^Mpc"^. Inserting this and M = 0.26Mq into equations ^ - and defining /_3 = //10"^Hz gives 



9iu(wd-wd) 



/ic(wd-wd) 



3 X 10^"/%^ , 
1.4 X 10~^Vr3''^ , 

4 X 10-^"/lJ/^Hz-^/2 



(18) 



These are respectively factors of 4, 2 and 2 less than found by the very detailed calculations of Schneider et al (2001). The 



differe nce is due largely to the different choice of normalization ar id is representative of the true uncertainties (see Phinney 
(2001) for an extensive discussion): Schneider et al used the Scalo (1986) IMF, a constant Galactic supernova rate of O.Oly"^, 
and a binary fraction of 100 percent. If we ad opted the same binary fraction, and scaled their Galactic birthrate (0.044y~^) 
to that of Iben, Tutukov & Yungelson ( 1997 ) (0.024y~^), the numbers would be in perfect (but spurious) agreement. The 
slight inaccuracies introduced by the theorem's neglect of systems whose merger occurs over a range of redshift (and our 
neglect of sources which exit common envelope evolution at frequencies much higher than 10~'*Hz) are small compared to the 
astronomical uncertainties. 

As a second example, we consider figure 4 of Rajagopal & Romani (1995). They computed the gravitational wave 



background due to mergers of super-massive black holes in galactic nuclei, at the wave frequencies accessible through pulsar 
timing / ^ 0.1 — ly~^. Their assumptions are approximately that 0.2 of bright galaxies contained a black hole, and that 
each bright galaxy had 5 mergers in its lifetime, or on average one merger with another black hole. Their black hole mass 
distribution has a black hole space density A'^o ~ 10~'*Mpc~"^ and typical black hole mass of ~ 10''-*Mq. Taking M for 
an equal mass pair, equation |l^ gives he ^ 8x 10~"(//ly~^)~^''^ However their mass function has A'o oc M up to 
3 X IO'Mq, and those heavier black holes merging with their typical ~ 10^ *Mq ones have M oc M" *, so the larger masses 
make contributions to oc M^'^^dlnM. Integrating over the mass function gives a total hc{f) ~ 2 x 10~^^(//ly~^)~^'^'', in 
very good agreement with the /si = 1 curve in the simulation shown in their figure 4. Of course, beliefs and data about the 
numbers and mass distributio n of su per-massive black holes have evolved considerably since 1995, so this value is not the one 
preferred today (see Phinney (2001)). 

It should be emphasized that in this discussion, we consider only the typical amplitude and spectrum of the backgrounds. 
The important issue of whether particular experiments could resolve the various backgrounds (by identifying and removing 
sources through optimal filtering of template waveforms, or by angular resolution, or both) is deferred to later papers. 



4 THE THEOREM AGAIN: MATHEMATICAL DERIVATION 

The flux of gravitational wave energy from a distant source is 

Denote the Fourier transforms of h+{t) and hx{t) by superscript tildes: 

/oc p OO 

/t+,x (t)e"""-'^' dt with inverse /i+,x (t) = / /i+,x (/)e'^"^'* dt . (20) 
-OO J —OO 

Apply Parseval's relation 

/OO 
df (21) 
- OO 

to the time integral of equation |l^, using the fact that the Fourier transforms of /i+.x equal z27r//i+.x (/): 
S{t)dt = / (hl+hl) dt 



2 G 



V- I f[\h+{f)\' + \hAf)\']df, (22) 



Theorem on gravitational wave backgrounds 5 



where in equation [22| we have made use of the fact that /i+,x are real, /i+,x(/) — ^+,x(— /), and hence \h+^x{f)\'^ = 
\h+.x (— ,f)P- This allows us to fold the negative frequency part of the integral onto the positive frequency part. If we average 
over source orientations f2s, or equivalently, over observer positions around a given source, 



(23) 



where Lg^ is the gravitational wave luminosity measured in the cosmic rest frame of the source, and d_L is the luminosity 
distance to the source. Since time t at redshift 2 = is related to time tr in the source's cosmic rest frame by dt= [1 -\- z)dtr 
we have 



47rdi 



where Eg^ is the rest-frame energ 

and 

1 



Comparing equations 

{S{t))n, dt 



1+2 

jitr) dtr — ~ ^o~Equ' , 

i-jidl 

emitted in gravitational waves. 
4, and using d/ = (1 + z)~^dfr we get 



(24) 



TT C 



2 G 1 

1 + z 
4ndl 



f{\h+if)\' + \hx{.f)f)^ dfr 



dEo 



dfr 



and thus identify 

'-di,f {\h+{f)\ 



dfr 



G 



(25) 



(26) 



where we have introduced the proper-motion distance dM ~ dL/{l + z) (cf. section 5 of Hogg (2000)), which is also l/(27r) 
times the proper ('comoving') circumference of the sphere about the source which passes through the earth today. 

Now consider sources undergoing the catastrophic events at redshift z, at rate A'' per comoving volume per unit of cosmic 
time tr local to the event. As seen from earth, in earth time dt, the number of events which occur in dt between redshift z 
and 2 -I- d2 is 

rf# _ 1 rfVc 
dz ' 



dt dz 



1 + 2 



(27) 



where the comoving volume element is (cf. Hogg (200C) and references therein) 

dVc . c 1 
dz 



(28) 



where E[z) was defined in equation The number of events which occur in a comoving volume between the cosmic times 
tr{z) and tr{z + d2) is 



Niz) = N— = iV i 

^> dz {l + z)HoE{z)- 

Thus equation ^ can be rewritten 

The energy density in gravitational waves at 2 = is then, inserting equation 

d# 



^g^{f)PcC df/f = 



{S{t))n, dt 



dt dz 



dz 



(29) 



(30) 



(31) 



Aixd^c 



dEqw 

fr 



dfr 



N{z) 



1 + 2 



fr- 



dEgu 

dfr 



N{z)cA-KdMdz 
df 



dz 



f 



which reproduces equation |j| and hence the statement of the theorem, equation 



(32) 



Notice that equation ^ and the result are independent of whether the timescale over which an individual source emits 
gravitational waves at the frequencies of interest is long or short compared to the observing time. The case when it is long 
compared to the observing time is appropriate e.g. for the two examples treated in section ^: double degenerate binaries in 
the LISA frequency band, lifetimes 10^-lO^y, and the super-massive binaries contributing to the pulsar timing background, 
lifetimes 10*-10V- In this case the number of sources contributing to the background is much larger than the number of 
events (binary mergers in these cases) in the observing time T, by the ratio Egjn/{LgwT). But the energy which each radiates 
in T is only a fraction (Lg^T)/ Eg^ of the total given by equation So the product of the number of sources and the energy 
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(spectrum) radiated by each is independent of T and I/gm. This must be so, since £gw is independent of whether the observing 
time is ly or lO^y. 

The observational prospects for removing the background by fitting individual sources do depend on the observing time, 
however, so it is instructive to re-derive equations ^ and |l^ for the case of merging binaries in circular orbits, keeping track 
of the number of sources contributing in each frequency interval, while explicitly showing the cancellation discussed in the 
previous paragraph. 

For circular binaries of separation a, masses Mi and A'h and quadrupole gravitational wave frequency = 2/Pt given 

by 

G{Mi + Afa) = n^fra , (33) 



and the energy flux at earth averaged over source orientations and orbital phases is 

■K ^ 2 ( AGHhK'h V ( fr V _ L,^ , . 

AG^ \ c^adM J \l + zj {l + z)2dlj' ^ ' 

If the rate of mergers per comoving volume is A^, then the space density N{a)da of systems with separations between a and 
a + da is given by the continuity equation 

^ (aiV(a)) = -NS{a) , (35) 

where we have assumed that most of the sources are born at much larger separations than those of interest. Then N{a) = 
N/{-a). Since 

d GM1M2 

= Lgw , (36) 



dt 2a 
one finds 



2a^L„^ , , -GMiM. 



GM1M2 



and N{a) = N^:^^ . (37) 



The number density of sources N{fr)dfr radiating at frequencies between fr and fr + df,. is related to N{a) by —N{a)da — 
N{fr)dfr, which gives, after manipulation 

^^■^"■^ " JgmW^ ■ ^^^^ 

The contribution to the specific intensity / (erg cm~^s~^sr~^Hz~^) of the gravitational waves from sources between z and 
z + dz is 

''-'■^^f^^^zW''- ^''^ 
Integrating over all redshifts gives the total specific intensity. Into this equation, insert the right side of equation ^ for S. For 
N(fr) use N{a)da/dfr, with N{a) from equation ^ for dVc/dzdO, use l/(47r) times equation ^ and recall that dfr/df = 1 + z. 
After using also equation [l3l one gets 



1 /.vrfi\ GM1M2 -da , 



(40) 



The quantity in parenthesis is N{z). The energy density per log (observed) frequency is, using / = /r/(l + z) in the second 
line, 



^an.if)PcC =^ Idn = —fl 
c I c 



N{z) GM1M2 (-d\na 



2a y d In /, 



2 f°° N{z) GM1M2 



3 L 1 + z 2 



dz 

2 1 1/3 



G(Mi + M2) 



7r2/3 (GX)5/3 ^2/3 p N{z)dz 



3 G f J, (l + z)V3- (41) 

The final form reproduces equation while equation ^ allows one to check the instantaneous number of sources contributing 
to a given frequency interval of the background. 
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5 CONCLUSION 

We have shown that in order to estimate the contribution to the gravitational wave background from a population of sources, 
one need not specify a cosmology. The only work involved is some relativity (often simple): computing frdEgw/dfr for the 
life-history of an individual source, and some astronomy: estimating the number of sources which have ever lived and died in 
a unit comoving volume of the universe, and some idea of at what redshift they did so. Insert the results in equation ^ to get 
the full spectrum of the background. 
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